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Abstract 

We develop algorithms for the numerical computation of the quadratic hedging strategy 
in incomplete markets modeled by pure jump Markov process. Using the Hamilton- Jacobi- 
Bellman approach, the value function of the quadratic hedging problem can be related to 
a triangular system of parabolic partial integro-differential equations (PIDE), which can be 
shown to possess unique smooth solutions in our setting. The first equation is non-linear, but 
does not depend on the pay-off of the option to hedge (the pure investment problem), while 
the other two equations are linear. We propose convergent finite difference schemes for the 
numerical solution of these PIDEs and illustrate our results with an application to electricity 
markets, where time-inhomogeneous pure jump Markov processes appear in a natural manner. 

Key words: Quadratic hedging, electricity markets, Markov jump processes, Partial integro- 
differential equation, Hamilton-Jacobi-Bellman equation, Holder spaces, discretization schemes for 
PIDE. 



1 Introduction 

In an incomplete market setting, where exact replication of contingent claims is not possible, 
quadratic hedging is the most common approach, among both academics and practitioners. This 
method consists in minimizing the L 2 distance between the hedging portfolio and the claim. Its 
popularity is due to the fact that the strategy is linear with respect to the claim, and is relatively 
easy to compute in a variety of settings. 

In its most general form, the quadratic hedging problem can be formulated as follows. Consider 
a random variable H E L 2 (.Ft,P) (which stands for the option one wants to hedge) and a set X 
of admissible strategies, which is a subset of the set of adapted processes with caglad paths. The 
quadratic hedging problem becomes 



minimize E 



tdS t - H 



over x S K and 9 € X 



(1.1) 



where S is a semimartingale modeling the stock price. If (x*, 9*) is a minimizer, we call 9* the opti- 
mal mean- variance hedging strategy and x* its price. This problem has been extensively studied in 
the literature, starting with the seminal works of Follmer and Sondermann ( 1986 ) andlFdllmer and 



Schweizer ( 1991 ) and until the complete theoretical solution in the general semimartingale setting 
was given in Cerny and Kallsen (2007). The case when S is a P-square integrable martingale is 



particularly simple and can be solved using the well known Galtchouk-Kunita -Watanabe decom 
position. The general case is much more involved, and has only been solved in Cerny and Kallsen 
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(20071 by means of introducing a specific non martingale change of measure (the opportunity 
neutral measure). 

The problem of numerical computation of the hedging strategy is an important issue in its 
own right, since various objects appearing in the theoretical solution (opportunity neutral mea- 
sure, Galtchouk-Kunita-Watanabe decomposition, Follmer-Schweizer decomposition) are often not 
known in explicit form. When the underlying asset is modeled by a Levy process, a complete semi- 



explicit solution was obtained in Hubalek et al. (2006) using Fourier methods. Their approach was 



extended to additive processes in Goutte et al. (|2011|). Laurent and Pham (|1999|) and [Heath et al. 



(2001 1 characterize the optimal strategy vi a an H JB e quation in continuous Markovia n stochastic 



volatility models while ICerny and Kallsen (2008) and 



Kallsen and Vicrthauer 



(20091 treat affine 



stochastic volatility models using Fourier methods. 

In this paper, we propose algorithms for the numerical computation of the quadratic hedging 
strategy in general Markovian models with jumps. We first review the HJB characterization of the 



value function, obtained in De Franco (2012). We only give a brief review, referring the readers 
to De Franco (2012) for full details and proofs because in this paper, we are interested in the 



numerical schemes for the computation of the hedging strategies and in applications to electricity 
markets. The value function of the quadratic hedging problem can be related to a triangular 
system of parabolic PIDEs, which can be shown to possess unique smooth solutions in our setting. 
The first equation is a non-linear PIDE of HJB type, but does not depend on the pay-off of the 
option to hedge (the pure investment problem), while the other two equations are linear. We next 
propose two finite difference schemes for the numerical solution of the linear and the nonlinear 
PIDEs. The convergence of these schemes is carefully analyzed and we provide an estimate of the 
global approximation error as function of various truncation and discretization parameters. For 
the numerical schemes, we concentrate on the infinite variation case, which is more relevant in 
applications. 

Our main motivation comes from hedging problems in electricity markets. These markets are 
structurally incomplete and often illiquid, owing to a relatively small number of market participants 
and the particular nature of electricity, which is a non-storable commodity. As pointed out in 



Geman and Roncoroni (20061 and Meyer-Brandis and Tankov (2008), due to these features, the 



electricity prices exhibit highly non-Gaussian behavior with jumps and spikes (upward movements 
followed by quick return to the initial level) and several authors have therefore suggested to model 



electricity prices by pure jump processes (Benth et al. 2007 Deng and Jiang 2005). 

On the other hand, since the spot electricity is non storable, the main hedging instruments in 
electricity markets are futures. A typical future contract with maturity T and duration d guarantees 
to its holder continuous delivery of electricity during the period [T, T+d]. Maturities, durations and 
amounts of electricity are standardized for listed contracts. This continuous delivery feature implies 
that even if the spot electricity follows a simple model, such as the exponential of a Levy process, the 
price of the future contract will be a general Markov process with jumps, non-homogeneous in time 
and spacep] T herefore, Fourier methods such as the ones developed in Hubalek et al. (2006) and 
Goutte et al.| ( |20lT ) cannot be applied in this setting. For this reason, in Section|5| after introducing 
a model for the futures prices, where the spot price is described by the exponential of the tempered 
stable (CGMY) or Normal Inverse Gaussian (NIG) Levy process, we derive the associated HJB 
equations and use the finite difference schemes to compute the hedging strategies and analyze their 
behavior. The numerical results illustrate the performance of our method and show in particular 
that the computation of the hedging strategies under the true historical probability (as opposed 
to the martingale probability, which does not require solving non-linear HJB equations) leads to a 
considerable improvement in the efficiency of the hedge. 

The paper is structured as follows. After introducing the model and the quadratic hedging 
problem in Section [2j we review the HJB characterization of the solution and the regularity results 
in Section [3] The finite difference schemes for the solution of the HJB equations, which are the 
main results of this paper, are presented in Section [4] In Section [5] these results are applied to a 
concrete hedging problem in electricity markets. The proofs of the convergence results are given in 
the Appendix A. 



If a single delivery length is fixed, it is possible to model the future price directly as the exponential of a process 
with independent increments, as in |Goutte et al.j ( |2011} . However this approach does not allow to treat problems 
involving futures of different durations, say, hedging a product with specified duration with the listed contracts. 
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2 The model and the quadratic hedging problem 

Let J be a Poisson random measure on [0, +00) x M defined on a filtered probability space 
(Q, J 7 , ^t,P), Tt being the natural nitration of J. We suppose that contains the null sets and 
also T = Ft where T > is given. Let also dt x v (dy) be the intensity measure of J where v 
satisfies the standard integrability condition J R (l A |y| 2 ) v (dy) < 00. We denote 

J(dt x dy) := J(dt x dy) - dt x v(dy) 

the compensated martingale jump measure. On this probability space we introduce a family of 
real-valued Markov pure jump process as the solution of the following: 

dZ]: 7 - :=u (r, Zl' z ) dr + [ 7 (r, z£f , y) J (dr x dy) , Z l t ' z = z,r>t (2.1) 
Jr 

for t G [0, T) and z G R. The asset price process S is given by S*' 2 = exp(Z^ z ). We make the 
following assumptions: 

Assumption 2.1. 
/C/- The coefficients. 

i) . There exists ~p > suc/i i/iai H/iH^ < ju. 

ii) . For all t e [0, T] and y GM. the functions z ^ u(t,z) and z — >• r y(t,z,y) belong to C 1 (R,R) . 

m). TTiere exist JQj > ; > and a positive locally bounded function p : R — > M + smc/i i/iat 
/or all y and all t G [0, T] we feare 

K*,z)-//(t,z')| < isr, c ip |z - z'| 

/X/- T7ie Levy measure. The Levy measures v{dy) verifies v(dy) = v(y)dy where v(y) := 
9{v)\y\ ~( 1+a ^ for some a G (1,2), where g is a measurable function bounded in a neighborhood 
of zero: 

0<m g < g(y) < M g , Vy G {-yo,Vo) 
for some positive constants m gi M g and some yo > 0. We also assume that 

lim g(y) = g(0~) and lim g(y) = .a(0+) with c(0 + ),.g(CT) > 0. 

y^0~ y->0+ 

Example: The tempered stable ( CGMY) processes, whose Levy measure is of the form 

\y\ 1+ L y<o + \ y \i+<* y>0 

for c_ > 0, c + > 0, A_ > and A + > satisfies the above assumption. 
[I]- Integrability conditions. The function 

r{y) :=max (sup (| 7 (*,*,!/) \,e^ z ^ - l) , p(y)^j 



verifies 



sup ^ <M and r G L 4 ({|y| > y }, f(dy)) 



o<|y|<yo \y\ 
[ND]- No degeneracy. The function 

T(y) := inf (Vr(*.*.w) - l) 2 ven/ses |r| := / r 2 (y)^(dy) > 

'• 2 v ' Jr 

[RG]- Regularity of the 7 function. 
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i). For any t,Z the mapping y — > y(t, Z,y) is twice continuously differentiable in a neighborhood 
of zero and there exist positive constants mi,m2 such that 

< mi < inf \^f y (t,z,y) \ and sup \j yy (t, z,y)\ < m 2 

t,z,\y\<Vo t,z,\y\<y 

In particular 7 as function of y is invertible on {~yo,yo): we call 7~ 1 (i, z,y) its inverse, 
ti). For all t,z G [0, T) x R, j y (t,z,0) = 1 

Hi). The function j y is Lipschitz continuous in the variable z: 

sup \j y (t, z + h,y)~ j y (t, z,y)\ < m 2 \h\ 

t,z,\y\<ya 

We denote K max := m&x(Kf lp , Kf ip ), 

ft :=pi + /V - 1 - i)v{dy) and ||/2|| := sup \fi{t, z)\ (2.2) 

Jr t,z 
In the rest of the paper we denote Hr]^ v := Jm>i T {y) v (dy) whereas ||t|| 2 v :— J R T 2 (y)iy(dy). 



It is well known that the SDE (2.1 1 has a unique strong solution ( Jacod and Shiryaev, 2003). 
On the Assumption [RG] — ii) 

Among the assumptions listed above, undoubtedly [RG] — ii) seems to be the most restrictive one: 
if for example the jump function is of the form 7(^,2, y) = j(t,z)y then the only possible choice 
would be j(t, z) = 1 for all t, z. In this paragraph we prove that this assumption could be relaxed 
by making a special change of variable. More precisely, given a process Z which does not satisfy 
this assumption, we look for a process L defined by L t = 4>(t, Z t ) for some smooth function (f> such 
that 

dip :=fj, L (s,L t /)ds + J 1 L (s,L t /_,y)J(dyd S ) (2.3) 

where /i L and j L satisfy Assumptions 2.1 in particular, d y j L (t, I, 0) = 1 for all t, I. If for such </>, 
the function z — > <f>(t, 2) is invertible and smooth enough to apply Ito's formula then 

7 L (i, I, y) :=Ht, </>-% I) + 7 (f, 4>-\t, 0,1/)) - I (2.4) 

n L (t, i) :=^(t, rHt, i)) + fi(t, 4r\t, i)At, r\t, i)) 

at oz 



\y\< 



In particular one has 
If we select for example 



i {l L {t^,y)- 1 {t,r 1 {t t l) t y)^{tA-\t,l))\v{dy) (2.5) 
ly(t,l,0) = ^(t ! r 1 (t,/)))7,(i,r 1 (*,0 ! 0) 



W,z):=[" ( f S m (2.6) 
Jo ly(t,s,0) 



then trivially j y (t,l,0) = 1 for all t, I. The following Lemma shows that this choice guarantees 
that the coefficients fi L and 7 verify Assumptions 



2.1 



Lemma 2.2. Assume that there exist positive constants mi,m.2 such that 

i). For all t, z 6 [0, T]xR the mapping y — » 7(i, 2, y) is differentiable at y = and 
< mi < \ j y (t, 2, 0)| < m 2 for all t,z€ [0, T] x K 
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ii). The function (t,z) — > r ) y {t 1 z, 0) is differentiable and 



dt 



Jy(t,Z,0) 



dz 



lv(t,z,0) 



< ni2 for all t, Z € [0, T] 



Hi). The function z — > 4f'Jy(t,Z,0) is Lipschitz continuous: 



|7,(M,0)-| 7y (M',0) 



< m 2 \z - z'\ for all t £ [Q,T], z,z' £ 



Then the functions fi L and 7 defined in (2.4|-(2.5| with the choice of <f> given by (2.6) verify 
Assumptions \2.1\ 



The proof of this Lemma can be found in De Franco (2012), Lemma 7.16 



The message of this Lemma is that, up to some regularity of the function 7 y at y = 0, it is possible 
to remove the assumption [RG] — ii): we could work with the process L instead of Z and derive 
all the results for L. By applying the function <f> we could then obtain the corresponding results 
for the process Z. We refer to Chapter 7 in De Franco ( 2012[ ) for further details. 

Nevertheless, we prefer here to work with assumption [RG] — ii) because it will make all 
computations easier to handle. 



Admissible strategies and the value functions 

To describe the set of admissible strat egies in the quadratic hedging problem we follow the ideas 
developed in Cerny and Kallsen ( |2007[ ): we first introduce the sets of simple strategies: 



V := jtf := X!>il[ ei , Ci+1 ), Yt £ L°°(J- ft ), ft < stopping times j 

v t := {ei [m I £ V) 

The set of admissible strategies is a subset of the L 2 (P) -closure of T>: 

X := Id £ V : [ 9dS £ L 2 (I 



We define the wealth process for all t, z, x by 



dXp z ' x ' e := 6 r -dSp z , X*' z ' x '' 



(2.7) 



(2. 



(2.9) 



where 9 represents the number of shares in the portfolio at time t. The set of admissible controls 
is then given by 



X(t,z,x) := jfl^T]! Be X x + J^ 9 r -dSp z £ L 2 (P)| 



(2.10) 



Consider a European option of the form / (Zt) where / is, for the moment, a bounded and 
measurable function. The quadratic hedging problem can be formulated as follows: 



(f(z^)-x^ e ) 



QH : minimize E p 

over 9 £ X(0, z, x) 
The dynamic version of QH gives us the value function of the problem: 



v f (t,z,x):= inf E p 
v f (T,z,x) = tf{z) - xf 



(/( 



ryt,Z\ ■y't,Z,X,i 

Zjrp j JY rp 



(2.11) 
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As remarked by several authors, the function Vf has the following structure 

Vf(t, z, x) — a(t, z)x 2 + b(t, z)x + c(t, z) 
In particular, taking / = 0, one has 



Vq ft, z, x) := x 2 inf E 

9eX(t,z,x) 



because the set X(t,z, x) is a cone. Consequently 



a (t, z) := inf E 



9 r _dSi 



(2-12) 



(2-13) 



9, — dS~ 



(2.14) 



This problem is known in the literature as the pure investment problem. The dual formulati on of 
this problem relates the function a to the so called variance optimal martingale measure (Cerny 



and Kallsen 2007 1. We recall here some fundamental properties on the function a, whose proof 
can be found in Chapter 5 of De Franco] ( |2012 ) . 



Theorem 2.3. Under Assumptions \2. 1[ [C , I, ND] the function a verifies 



- C{T - t] <a(t,z)<l 



where 



C 



Furthermore, there exists T* > and Kf ip > such that if T < T* then 

\a{t,z')-a{t,z)\<K? ip \z-z'\ 



for all t G [0, T] and z,z' G K. T* depends on fj,, t, C and K max defined in Assumptions 2.1 
Moreover T* — > +oo when K max — > and the other constants remain fixed. 

Remark that these results hold true without assuming any particular structure of the Levy 
measure v(dy). The next goal is to characterize the functions a, b and c as the solutions of certain 
PIDEs. 

3 HJB formulation and main regularity results 

Remarks on notation For a function / : [0, T] xR —> R we denote WfW^ :— sup t<Ta , eK \f(t,x)\. 
For a function ip defined on [0, T] x K and k € N wc denote D k ip :— d k tp/dx k whereas dt<p denotes 
the derivative in the time variable. We adopt the following convention: for any I £ K + 

l=[l\+ {l}~, where {I}- £ [0, 1) 
I = \l] +{?}+, where {/}+ G (0,1] 



Let us first introduce the functional spaces in which we will work: for j3 G (0, 1] we define 

(W P > ■= sup ^ , (ip)fa := sup 



((3) _ \Tp{x + h) -ip(x)\ (f}) _ \<p{t,x + h) - ip(t,x)\ 

.— sup , vPIQt " — P \u\f 

xfl<\h\<l \n\ p * t,x,0<\h\<l l«r 



1/3 



The elliptic Holder space of order /, H (M. n ), is defined as the space of continuously differentiable 
functions tp for all order j < |~f| with finite norm 



i^iii : = EE \\^\L + E< D ^) m) 

j=o (j) (m) 



(3.1) 



G 



where Y^m represents the summation over all possible derivative of order j. The parabolic Holder 
space H l ([0, T] x R") is defined as the set of measurable functions <p : [0, T] —> H l (W l ) with finite 
norm 

m 

imi« : = EE W^ML + E<^C ( 3 - 2 ) 

j=° y) cm) 

The spaces defined above are all Banach spaces equipped with their respective norms. For a com- 



plete description see for example Chapter I in Adams and Fournier (2009). 



In the spirit of HJB approach we now introduce the operators associated to the process Z: 

Definition 3.1. For a real valued function ip e H a+s ([0,T] X R), 6 > 0, we define the following 
linear operators 

A t <p{t,z) := - (i—(t,z) 

Bt<p(t,z) := j \ y{t,z + i(t,z,y)) - <p(t,z) -j(t,z,y)^(t,z)J v{dy) 
Qt<f(t, z ) z ) + / ( e7 - 1) (<P(t, z + l(t, z, y)) ~ <f(t, z)) v{dy) 



Gt<f{t, z) := J (e 7 - l) 2 tp(t, z + j(t, z, y))v{dy) 

where fi stands for fi(t, z) and so on. In addition, % denotes the nonlinear operator 

Ht[tp](z):= inf [2irQ t tp(t, z) + n 2 G t <p(t, z)} 

|7r|<n 



whe 



CT 

n :=^max (HAlL^^Mli,, + ||r||y) (l + K? ip ). (3.3) 



|r, 

The main result concerning the functions a is 



Theorem 3.2. Let Assumptions \2.1\ hold true and consider T < T* as in Theorem \2.3\ The 
function a is the unique solution of 

da 

= -Q t +A t a-B t a-Ht[a}, a(T,z) = l (3.4) 

in the Holder space H a+S ([0, T] xR) for < 8 < a— 1. The function 1 1— > a(t, z) is also differentiate 
on (0,T). The optimal strategy for the stochastic control problem (2.131 is 

6* = e - Zt -ir* (t, Z t -) Xf_, Xf :=x+ [ 0* r _dS r 

where 







Qa{t,z) 
Qa[t, z) 



□ 



For the general value function Vf we have 
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Theorem 3.3. Let T < T* as in Theorem \2.3\ Let also Assumptions \2.1\ hold true and f G 



H a+d (R) for some < S < a — 1. The function Vf in ( |2.1l[ ) admits the decomposition 

vj(t, z, x) = a(t 7 z)x 2 + b(t, z)x + c(t 7 z) 
where a is defined in (2.14), so it does not depend on f , and it is the unique solution in H a+S ([0, T] x 



of 3.4 whereas b and c are the unique solutions of the following linear parabolic PIDEs 
db 



= - 




dt 

dc 



Ab- Bb-n*Qb, 
1 (Qb) 2 



+ Ac~Bc , 
at 4 ya 



b(T, .) = -2/; 
c(T, .) = f 2 



(3.6) 
(3.7) 



in the Holder space H a+s {[0,T]xR), where 7r* is defined in ( |3.5[ ). The functions 1 1— > a(t, .),b(t, .),c(t, .) 
are also differ entiable on (0,T). 

Furthermore the optimal policy in the control problem (2.11) is given by 

0* — c - z t- 



n*(t,Z t _)Xl- 



1 Qb(t,Z t -) 
2Ga(t,Z t _) 



XI :=x- 



01_dS r 



(3.8) 



□ 



sition (2.12) we also obtain the optimal price in (2.11): 



The proof of these results can be found in Chapter 7 of De Franco (|2012). From the decompo- 

(3.9) 



x (/ := arg wiv f {t,z,x) = - —7- 

xer 2a(t, z) 



which is a linear function of the payoff / since W is. 



Non smooth payoff Theorem 3.3 allows us to characterize the value function Vf when the payoff 
function / is sufficiently smooth, i.e. / e H a+S (R). However, in most cases of interest (for example 
put options, straddles or bear spreads) this function is not even continuously diffcrentiable. The 
following lemma proves the stability of the optimal price x*(f) and the optimal hedging strategy 
under small perturbations of the function /: 

Lemma 3.4. Let f\,fi be two measurable functions with fi(Z!f z ) G L 2 (P) for all t,z,i = 1,2. 
Then for any t < T and z G M 

\x*(f 1 )(t 7 z)-x*(f 2 )(t,z)\ <a(i, z)-^2 I) (/l _ /2 )(4*)|| 



IK 1 - v^) (t,z,x)\ <2 (x + || (A + / 2 )(4 z )|| L2(p) ) || (/1 - j2)(4 Z )L 



(P) 



Fix now (t,z,x) and let f n such that ||(/ n — f)(Z^ z )\ 



> 0, n — > 00. // 9 n is the optimal 
control in the problem (|2.11|) with payoff function f n then, for all e > 0, there exists some N > 



such that for any n> N one has 

v f {t,z,x)-E F 



[f{Z l T z )-x- J" 



0™ dS~ 



< e 



The proof of this Lemma can be found in Chapter 5 of De Franco (2012). One can thus 



approximate a non-smooth payoff function / with smooth functions /„, controlling the error on the 
value function and the cost of the hedging strategy with ||/ — /„|| 2 . Furthermore the corresponding 
strategies (8 n ) become e-optimal for the pay-off / starting from sufficiently large n. 
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4 Numerical solution schemes 



We now present a numerical scheme to solve the PIDE introduced in Section [3] when the Levy 
measure v verifies the Assumption 2.1 [L]. From (3.8) and (3.9) we remark that in order to solve 



the problem (2.11), i.e. to find the optimal strategy 9* and the optimal price x* , we only need to 



compute the functions a and b, solutions, respectively, of PIDEs (3.4) and (3.6) 



In order to solve PIDEs ( 3.4 ) and ( 3.6 ) numerically, we construct a finite difference discretization 



scheme using the Markov chain approximation technique developped in Kushner (1976). One of 



the advantages of the probabilistic treatment is that it allows to estimate the error due to the 
truncation of the domain in a simple manner. 

We calculate the values of a and b on a regular space grid Zj = jAz for j 6 (— N, N) and 
regular time grid ti = iAt, i = 0, . . . , Nt with At = The values of the approximations of the 
functions a and b at points (ti,Zj) will be denoted, respectively, by a* = a l {zj) and bj = b l {zj). 
For j (£ (-N, N), the approximation will be defined by the boundary condition: a* = q a (ti 



and 
and 

q b (T, Zj) for all j will stand for the final boundary condition. In the numerical examples, we take 
q a = 1 — this choice of boundary condition can be justified by the fact that if S is a martingale, this 



bl = q b (ti, zj), where the functions q a and q b are assumed measurable and bounded. q a (T, Zj 



is indeed the exact solution. In Theorems 4.2 and 4.8 we shall see that the effect of the boundary 



conditions on the approximation becomes negligible for AzN sufficiently large. 

We now detail the computation of the integral terms of the intcgro-differcntial operators in 



(3.4). In order to avoid interpolating the values of a and b, we use a space and time dependent grid 
for discrctizing the Levy density, with discretization points denoted by (yi(t, z))—m<i<m- We select 
the discretization point yi(t, z), which corresponds to the center of z-th discretization interval, as the 
unique solution of the equation j(t, z, yi(t, z)) — iAz. The boundaries of the discretization intervals 
will then correspond to half-integer values of i. Although these discretization points depend on t 
and z, we will sometimes omit this dependence to simplify notation. 

To treat the singularity of the Levy density at zero, we adapt the methodology inspired by 



Forsyth et al. (2007) and divide the real line into four disjoint regions, for an integer k > 1 



^o(M) = [y ■ V- K -i 2 (t,z) < y < y K+ i(t,z)} , 

toi(t,z) = jy : y K+ i(t,z) < y < y c+ i(<,z)| U \y : y_ K _i(t,z) >y> y_ c _i(t,z)} , 

Sh{t,z) = [y ■ <y - 2/m+a(M)} U jy : y_ c _i(t,z) >y> y_ M _±(t,z)j , 

toa(t,z) = [y ■ V < I/-m-a(M)} U |y : y > y M+ i(t,z)j , 

where £ = inf{z : j/j(t, z) > 1} and Q' = inf{i : y-i(t, z) < —1}. Without loss of generality, we shall 
always assume that k < £ < M and n < £' < M. 

The jumps of the Levy measure in the regions Oo, . . . , ^3 are treated as follows: 

• The small jumps in the region f^o are truncated and the corresponding part of the integral 
operator is replaced by a local operator (e.g., second derivative). 

• The jumps in the regions f2i and ^2 are discretized. 

• The large jumps in the region ^3 are truncated. 

For reader's convenience the different truncation and discretization parameters of our algorithm are 



listed in the following table. Theorem 4.2 shows how these different parameters affect the overall 
approximation error. 

Parameter Meaning 

At Time discretization step 

Az Space discretization step 

N Space grid size 

M Levy measure truncation (number of points) 

k Truncation of small jumps (number of points) 
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Set Ay,-(i, z) = y i , i (t, z) — y i— 1 (i, z) and introduce the weights : 



1 



2/i+l/2(M) 



y 2 v(y)dy if g/* € fii 



y>-i/2(M) 



v{y)dy 







if J/i G ^2 



otherwise 



(4.1) 



The integrals z) can often be calculated explicitly, such as for CGMY process, where they 
can be expressed in terms of the incomplete gamma function. Otherwise, they may be calculated 
numerically, using, e.g., the 5-point trapezoidal rule, which, in the second case above yields 

In the following, we shall assume that these integrals are calculated explicitly without error; the 
additional error introduced by the their numerical evaluation can be easily estimated along the 
lines of other computations in the following section. 



Definition of the approximating Markov chain Let the weights a and j3 be given by 

D(t, z) fi(t, z) 



a(t, z) 



2Az 2 



2Az 2 

if both these expressions are positive or and by 

D(t,z) 



2Az 
u(t,z) 
2Az 



(4.3) 
(4.4) 



a(t, z) 



2Az 



max(0, 



(i(t,z) 



n , s D(t,z) . u(t,z). 

0(t,z) = 1 ^+max(0, / ^i) 

otherwise, where the coefficients D and fi are defined by 



Az 

^z 
Az 



(i(t,z) = n(t,z) - ^2 cj l (t,z)7(i,z,y,(t,z)), 

K<\i\<M 

and D(t,z):= / j(t, z,y) 2 v(y)dy. 



(4.5) 



Assume that At is small enough so that for all (t, z 
D(t,z) fl 



^2 \i\ui(t,z)+ ^2 < 



Az 2 Az ' ^ ' ^ At 

K<\i\<M K<|i|<M 



(4.6) 



We construct a discrete-time process (■Zt j )i=i ) ...,jv T approximating the process ( |2.1| as follows: 

' Z ti + jAz,j = k + 1, . . . , M, with probabilities u)j(t, Z ti )At 
Z ti + jAz, j = —M, . . . , — k — 1, with probabilities Uj(t, Z ti )At 
Z u +Az with probability /3(t,Z u )At 
Z ti — Az with probability a(t,Z ti )At 

Z u with probability 1 - a(t, Z u ) + /3(t, Z u ) + u>j(t, Z u ) At 

\ K<\j\<M J 



Zt i+1 = < 
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When the chain starts from the point z at time t n , its value will be denoted by Z^ tn , i > n. 
To make the notation more compact, we shall denote the probability of transition from Z t . to 
Z u + jAz by Pj (t h Z u ), for —M < j < M. 



Remark 4.1 (On the CFL condition). The condition (4.6) ensures the positivity of the transition 
probabilities of the Markov chain and hence the stability of the resulting scheme. By definition of 
the weights D and u>, 

1)11 "'' '" Yl z2 



Az 2 Az 



1 



Az 2 
1 

Az 2 



K<\i\<M 

j(t,z,y) 2 v(dy) + 



l 



K <\i\<M 

l(t,z,y)\v(dy) 



v(dy) 



Oiuo 2 



JJ_ 
Az 



|u|<Az(K+i) 



ii 2 i/o7 l (du) + 



Az(K+h)<\u\<MAz 



Az 



+ 1 1/07 L (du) + 



Az 



Since under our assumptions j y (t,z,y) ~ 1 when y — > 0, this can be further approximated by 



1 

Az 2 



M<Az(K+i) 



u 2 v{du) 



(c+ + c^)Az" 



Az(k+A)<|u|<A/Az \Az 
\2-a 



+ 1 I/(dtt) + 



A;. 



a (( K +l/2) 2 ~ a (k + 1/2) 1 -" (k + 1/2)-" 



2- a 



a — 1 



A j 



Hence, the stability of the scheme requires At < 



4.1 A finite-difference scheme for the function a(t,z) 
Let 

= inf{i > n : Z?f n <£ [-NAz, NAz}}. 
Consider the following discrete-time control problem for Z: 



a n (z 3 ) = inf E 



P tn - Z i/\N T 

n (1+^^ 



1)) q «(f3^AN T ,Z^ ANT ] 



where 7Tj € [—11, II] for i = n + 1, . . . , 7V T . 

The dynamic programming principle for this discrete-time control problem writes 



a n { Zj ) = inf _ E 



i + 7r ( e ^ +1 -^-i) a n+1 (z;;;/) 



, if j € (-N,N), (4.7) 



7rG[-n,n] 

a"fe-) = 9°(^ J ^), if j<£(-N,N), 
or in other words 

0?= inf. £ m {t n ,z j ){l + -K{e lAz -l)faj;Xl 
7re[-n,n] _ M < ; < M 

which gives a fully explicit finite difference scheme for approximating the function a. Using the 
explicit form of the transition probabilities, it can be rewritten as 



At 



(a(t n , Zj ) + (3(t n , z 3 )) a] +1 + [3(t n , Zj)a!£l + a(t n , Zj )a^ + £ w,(*n, Zj)(c$$ - a] +1 ) 



K<\l\<M 



inf 

7re[-n,n] 



tt 2 ]T ^ Z J)^ AZ - l ) 2a ltl + Zj)(e Az - + a{t n , z^e^ - 1) 



2 n+1 



lYa 



3-1 



. K<\l\<M 



-Az 



l)a 



n+1 



lK<|i|<Atf 



(4.8) 
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Hence, this scheme uses the following approximations for the operators A, B, Q and Q appearing 
in Equation (3.4): 

(B t - A t )a(t n , zj) « - {a{t n ,zj + (3(t n , zj) a] +1 + /3(t n , zjtfjl + a{t n , zjajj\ + ^ u)i(t n , Zj 

K<\l\<M 

g t a(t n , Zj ) « J2 ^)( e ' Az - i) 2 ^ 1 + ^)(e Az - l) 2 a^i + a(*n, *i)(e" A " - l) 2 a"+i 

K<|i|<M 

e t o(t n , zj « ^ «i(t„, z,)(e iAz - l)^ 1 + /3(i„, z,)(e Az - l)^ 1 + ^(e^ 2 - \)d!*\ 

K<\l\<M 

The following theorem analyzes the error and convergence of this algorithm. 

Theorem 4.2. Assume that the following technical conditions hold in addition to Assumptions 
DO' 

1. The Levy density v is twice continuously differentiable with bounded derivatives outside any 
neighborhood of zero; 

2. The function g has bounded derivatives near zero; 



3. Assumption 2.1-\RG\ — i) holds true with yo = +oo; 

4- The function j(t,z,y) is 3 times continuously differentiable with respect to y with bounded 
derivatives; 

5. The data of the problem are such that the functions a and b together with their derivatives 
up to order 4 with respect to z and up to order 2 with respect to t are bounded. 



Let the condition (4.6| be satisfied and assume that NAz > MAz + \xT . Then, there exist two 
positive constants c, C < oo, which do not depend on truncation of discretization parameters, such 
that for all k > c, 



(i + M) 



^3 — a t ,3 — a , \ ^3— a. ^1 — a 



\a n (zj - o(tn, zj\ < C<j {N 2 M )Az-fiT +At + AZ K + A 

(l + \y\+T(y) + T 2 {y))v(dy) \, 

(| v |> 

for -N <j<NandO<n< N T ; 

The proof of this result is given in Appendix A. 

Remark 4.3. This theorem gives the decomposition of the global approximation error into four 
sources: the domain truncation error (first term in brackets), the time discretization error (sec- 
ond term), the error coming from space discretization and truncation of small jumps (third and 
fourth term) and the big jump truncation (last term). It may be used for the optimal choice of the 
parameters of our algorithm. 

Remark 4.4. The constant c is needed to ensure that for Az sufficiently small, the weights a 
and j3 may be assumed to be given by (4.3 1 -(4.4 1. If the Levy measure is locally symmetric near 
zero, that is, c+ = c_, one can take c — 1. Indeed, in this case, the integral J E<X<1 xv(dx) remains 

bounded as e —¥ 0, which means that /t(i, z) remains bounded as well. On the other hand, D ^jjji is 

of order of {nAz)~ a , which means that, eventually D ^jj 3 becomes bigger than for all n > 1 

(since under our assumptions a> I). 



12 



Remark 4.5 (Implicit-explicit scheme). In practice, it may be preferable to use an implicit scheme 
for the convection diffusion part and an explicit one for the integral part: 



At 



(a(t n , Zj) + f3(t n , zj)) a] + /3(t n , Zj)a] +1 + a{t n , Zj)tf_ x + ^ u>i(t n , z 3 ){a™+l - a^ +1 ) 



«<|2|<Af 



inf 

7re[-iLn] 



, K<|2|<M 



yK<\l\<M 



(4.9) 



The error analysis for such a scheme may be done by introducing a two-dimensional Markov chain 



step requires that at least 



along the lines of (Kushner 1916, section 6.5). We observe however, that the stability of the explicit 

J2 w ^' z )-at 



K<\i\<M 



which still leads to a condition on At of the form At < for some constant C . 

Computing the optimal strategy The optimal hedging (pure investment) strategy may be 
computed as the value of the maximizer in (4.7): with 



7T*" =7T* n (Zj) =-IIV 



V 



E 



E 



e '■n.+ l 



ail 



The error of evaluating the hedging strategy may be greater than that of the evaluating the value 
function. The following corollary gives an estimate of this error. 



Corollary 4.6. Under the assumptions of Theorem 4-2, 



< 



c (i + M) 

/Al\(N-M)Az-p,T 

(l + \y\+T(y)+T 2 {y))v(dy) 



At + Az 3 - a K 3 - a + Az^k 1 -" (4.10) 

(4.11) 



The proof of this Corollary can be found in Appendix A. 

Remark 4.7. For the error of approximating the optimal strategy to tend to zero, the space dis- 
cretization step Az must therefore be sufficiently small compared to At. At the same time the CFL 
condition imposes a lower bound on Az, which may not tend to zero faster than At 1 /". Letting 
Az ~ At 1 /", we get a convergence rate of At2 A (a~i) for the optimal strategy. 

4.2 A finite-difference scheme for the function b 

The finite-difference approximation for the function b is defined as the solution of the following 
linear dynamic programming problem: 



b n (z j )=E[b n+1 (Z t t :^)]+^ n (z j )E 



«»+* MC) 
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for j e (-N, N) and b n (zj) = q b (t n ,z j ) for j £ (—N,N). In other words, 

b l= E Pl(tn,Zj)(l+7r* n (e lA *-l))b£t. 

— M<KM 



(4.12) 



The terminal condition is given by bj T = —2f(zj) and in the numerical examples we take 
q h (t, z) = -2 /Q 7 ) as well. 

Equation (4.12 1 defines a fully explicit finite difference scheme for computing the function b, 
which can be rewritten as 



b l +1 b l 
At 



(a(t», Zj) + /3(t„, z 3 ))V] +1 + /3(t„, z^+l + a(t n , zj)a]+l + £ w,(t», - b] 



K<\l\<M 



+ 7T*" [ ^ Wl (*„, ^)(e /Az - l)^^ 1 + P{t n , Zj )(e Az ~ \)V*\ + a(t n , Zj )(e- A * - 

,K<|i|<M 



(4.13) 



This scheme uses the same approximations for the operators appearing in (3.6) as the scheme for 



the function a defined in section [4. 1| An implicit-explicit scheme for the function b can be defined 
along the lines of Remark |4.5| 

Since the numerical scheme for the function b uses the approximate pure investment strategy 
7r* , the error of approximating b is determined by the error of approximating the optimal strategy. 



Theorem 4.8. Under the assumptions of Theorem \4-2 
\b n ( Zj )-b(t n ,z 3 )\< ° J " 



At j (N - M)Az - p.T 

{l + \y\+T{y)+T 2 {y))v{dy) 



At + Az 3 ~ a K 3 - a + Az i - a K 1 ~ a 



The proof of this result is provided in Appendix A. 



5 Application to the electricity market 



Many studies have shown that price spikes in electricity and gas markets are incompatible with 
Gaussian dynamics (Geman and Roncoroni 2006 Meyer-Brandis and Tankov 20081 and and 



several models based on Levy processes have been developed to fit the observed fat tails (Deng 



and Jiang 


2005 


Benth et al. 


2007 



or make some assumption on the change of probability resulting in a similar model under the 



martingale measure ( Benth et al. 2007 1 . However, these martingale models are not adapted for 
the evaluation of hedging strategies, since the hedging error should be computed under the historical 
measure. 

In this paper, we propose a model which describes the deformation of the forward curve directly 



und er the historical probability, and satisfies Assumptions 2.1 so that we can apply Theorems |3.2 
3.3 We start by introducing a Levy process L as follows 



L s = £s - 



yj(ds x dy) 



(5.1) 



where ( S 1 and J is a compensated Poisson random measure, whose Levy measure is denoted by 
v(dy). Fix c e R+, l{s) = e~ cs and 



s dL s 



(5.2) 
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We model the price at time t of the future contract with maturity T and instantaneous delivery 
as a random perturbation of the initial forward curve tf> which is supposed to be known. Using the 
above notation we have 

F ,T,t = mT)e l(T ~> A t 

By no arbitrage in the futures market, the price at time t of a future contract with duration d is 
equal to the average over the time period [T, T + d] of the future contract prices with instantaneous 
delivery. We therefore model the price at time t of a future contract with delivery time T and 
duration d > by 

1 r T+d _ i r T+d 

F d>T>t =jJ t Fo^tds =-j V(0, s)e L{s)At d S 

For reasons which will become clear in the sequel, we prefer the following notation: 



F d ,T,t ■= exp($(A t )) where 9(A) := log (i J + V(0, s)e l ^ A ds J (5.3) 



The model (5.3) essentially states that the price of a future contract F dt x,t 1S the average price 
on the interval [T, T + d] of the future contract with instantaneous delivery up to the random 
perturbation e l ^ A . In this context, the problem of hedging a European option on Fd,T,t with the 
quadratic hedging approach becomes 



minimize E 



^H(F diT ,t) - x - jf 8 u -dF d , T ,^J 



over 9 and i£R (5.4) 



for a given map H. The process F di T,t corresponds to S in the formulation (1.1). The following 



results proves that Z = log(-F) is a Markov jump process satisfying our assumptions. 
Lemma 5.1. The process Z t := log(i 7 d,T,t) verifies 



dZ t = fi(t, Z t )dt + J j(t,Z t -,y)J(dydt) 

where 

j{t,z,y) :=$($- 1 (z)+ye c *)-z 

»(t,z) :=Ce ct $'(a>- 1 (z))+ f ( 1 {t 1 z,y)~ye ct <$>'(<S>- 1 {z))) V {dy) 



Assume that the Levy measure v(dy) is given by v(dy) = g(y)\y\ ^ 1+a ^dy, for some a £ (1,2) and 
a bounded, strictly positive and measurable g such that the following conditions hold true: 

i) . There exists m < oo such that for all y, y' £ (~yo, 0) U (0, yo) with yy' > 0, \g(y) — g(y')\ < 

m\y - y'\ 

ii) . lim g(y) — g(0~) and lim g(y) = g(0 + ) with g(0 + ), g(0~) > 

y->-0- y-*0+ 

Hi). I y i v(dy) + / e Ay v(dy) < +oo 

Jy<-1 Jl<y 

Then the functions /i and 7 verify the Assumptions \2J\ [C , L, I, ND. RGi, RGiii] , where the func- 
tion t is given by 

r(y) = max(|y|, \e v - 1|) 
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The proof of this result is given in Appendix B. 



In order to apply our results (Theorems 3.2 and 3.3 1 we also need to verify Assumption 2.1 



[RGii] and it is easy to see that the function 7 does not verify it: however, as we have already 
said in Section [2] this can be avoided by making a change of variable L t = <f>(t, Z t ). We refer to 
Chapter 7 and 8 in |De Franco (2012) for further details. 



In terms of the process Z, the problem (|5.4|) becomes 
Vf(t,z,x) 



inf E 

6eX(t,z,x) 



(5.5) 



where X(t, z,x) is defined in (2.10) and f(z) — H(e z ). 



We now describe a special class of pay-offs which is of interest in problem (5.5). Let p(x) 
(G — x) + and define 

1 f T ' 



iP(0,s)e lis)A ds 



for d! ^ d. From ( |5.3[ ) it follows that h o $ _1 (Z i ) = Fd>,T,t, and then, by defining / := p(h o 

we obtain f{Z t ) = (G — Fd>.T,t) + > which is a put option written on a future contract with different 

duration d' . Using this specific option we can rewrite problem (5.5) as follows 



inf E l ' z ' x 
9eX(t,z,x) 



(G - F„/ iTit )+ - x - 



0,,-dF, 



d,T,u 



The financial meaning of the above problem is particularly interesting: one tries to hedge (in the 
quadratic sense) a put option written on a future contract with duration d' 7^ d using as hedging 
instrument the future contract with duration d. This may be useful when, for example, one sells 
a future contract with a non-standardized duration in the OTC market and hedges the resulting 
position using instruments which are liquidly traded. 

5.1 Numerical Example: the CGMY model 



In this section we study the problem (5.5) when L in (5.1 1 is a CGMY process (Carr et al. (2002)) 
We can write then 



L t = (fi + CT(1 - Y)(M Y - 1 - G Y - 1 )) t + 



where J is a compensated Poisson random measure with intensity 



yj(dyds) 



e -My 

v{dy) = v{y)dy, v{y) = C—^yl y>a + C 



aGy 



J^Y 1 y<0 



The goal of this paragraph is to solve numerically the equations (3.4) and (3.6) when the model 



for Z is given in Lemma 5.1 and the source of randomness L in (5.1) is given by the CGMY 



process introduced above. We first apply the implicit-explicit variant (4.9) of the scheme (4.8) for 



the function a with maturity T = 1. The coefficients /i and 7 given in Lemma |5.1 
fx(t, z) =§\<S>- l ( z )) ( M + CT(1 - Y)(M Y ~ 1 - G Y - 1 )) e ct 

+ [ U(^- 1 (z)+ye ct )-z-ye ct ^'(^- 1 [z))] v{dy) 
Jr l - 1 

l{t,z,y) =$ {^-\z) +ye ct ) - z, 



and the function D(t,z) introduced in (4.5), are computed numerically. To this end, we compute 
the integration points y^ as in Section |4J such that 7(t, z, yi(t, z)) — iAz, or equivalently 

Vl (t, z) := e- ct + ^Az) - ^(z)) 
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We choose k = in the definition of fio(i> z )i an d by expanding 7 around zero, we obtain 



D(t, z) := / 7 (i, z, y) 2 ^) =s e 2ct ($ (^(z))) 2 / y 2 ^) 

Using an approach similar to the one of Section [4j we approximate 

$ + ye ct ) - 2 - j/e rf $'($ _1 («))l u(dy) 

with 

iAz($'($- 1 (z))e ct ) 2 ^ + ^ d>,(t,«) [$ ($- 1 (^) + 2/«e c4 ) - z - 2/i e ct $'($- 1 ( 2 )) 

w,l«l>* 

= iAz($'($- 1 (z))e ct ) 2 ^ + ^ <&<(M) [iAz- (^- 1 (z + iAz)-^- 1 (z))^'(^- 1 (z)) 

Vi, \i\>k 



where the weights t2>i(i, z) are given in (4.2). The effect of approximating the coefficient of the 
PIDE on the solution of the source problem, appearing for example in the pricing of European 



options, has been studied in Jakobsen and Karlsen (2005) 



We solve the problem (5.5) for European options /, with maturity one week and delivery for 

itract in this cas 

i>(0,s)e l( - s)At ds 



the 7 days of the week. We recall that the future contract in this case is given by 

/" 14 

F7days,lweek,t 



1 

7 J 7 



(5.6) 



and At is given in (5.2) with L being the CGMY process defined above. The initial forward curve 
for the seven days of delivery is given in Table [l] indicating that prices are lower for the week end. 
This case with continuous long delivery corresponds to a non stationary process where the hedge 



Day 


s 


Price O(0, s)) 


Monday 


se [7,8) 


80 


Tuesday 


s e [8, 9) 


90 


Wednesday 


s e [9,10) 


70 


Thursday 


s e [10,11) 


90 


Friday 


s e [11,12) 


80 


Saturday 


s G [12,13) 


70 


Sunday 


s G [13, 14] 


60 



Table 1: The forward curve. Prices are given in Eur. 



cannot be calculated efficiently as in Hubalek et al. (2006) or Goutte et al. (2011) 



We use the scheme (4.9) to obtain a numerical approximation of the function a and the optimal 
control tt*. The trend is fi = 0.01, and the CGMY parameters are Y = 1.2, C = 0.01, G = 1.1, and 
M = 1.1. The space mesh step is always chosen such that N Az = 8, and the number of points M 
for approximating the Levy measure is equal to ^ . 

First we analyze the the convergence of the approximation as function of the number of time 
steps and the space mesh size. Table [2] shows that effect of increasing the number of time steps 
on the computed value is less pronounced than that of increasing the space mesh size, which is 



consistent with our main theoretical result (Theorem 4.2 which shows that the convergence as 
function of At is faster than that as function of Az). 

In the remaining numerical experiments in this section we suppose that N = 800, M = 200, 
take the mesh step equal to Az = 10~ 2 and the number of time steps equal to 800. Figure [TJ shows 
the function a whereas the optimal control it* is given in Figure [2] Finally, Figure [3] shows the 
(opposite of the) optimal price b when we consider an at-the- money put option in (|5.5[) . We observe 
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Time steps (N T ) 


Space steps (N) 


Price (moneyness 1) 


800 


800 


4.91 


400 


800 


4.88 


200 


800 


4.83 


100 


800 


4.73 


800 


400 


5.19 


800 


200 


5.63 


800 


100 


6.32 


800 


50 


7.26 



Table 2: Convergence of the optimal price of the quadratic hedging strategy as function of time / 
space discretization parameters for the CGMY model. 



that the value of the function a (cost of optimally hedging the pay-off equal to 1) increases as the 
strategy approaches maturity, but has little sensitivity to the price level of the underlying asset. 
The magnitude of the optimal strategy it descreases with time, which is explained by the fact that 
it becomes too risky to invest a lot when the trading horizon is close; minimizing the quadratic 
functional is equivalent to maximizing a concave utility — x 2 , which means that the investor is risk 
averse. The magnitude of the function b also decreases with time: options loose value when they 
approach maturity. 




Figure 1: The value function (t, z) i-> a(t, z) for the CGMY process. 

Practitioners usually price options of this type and calculate the hedging strategy assuming 
that the underlying process F is a martingale. It is therefore interesting to evaluate the loss of 
efficiency when using the hedging strategy computed in the martingale model. Assuming that F 
is a martingale means that we should have 

! rT+d 

F d ,T,t -=2] ^(0, s) exp (M (s, t) + l{s)A t ) ds 
for some M that makes F a martingale under the historical probability P. Using Lemma 15.1 in 
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Pi curve 




Figure 2: The optimal control (t, z) ir*(t, z) for the CGMY process. 



B function for put 




Figure 3: The value function b(t, z) for an at-the-money put option for the CGMY process. 



Cont and Tankov ( 2004 1 we obtain 

t 



M(t, s)=- (fi + CT(1 - Y)(M Y - 1 - G Y ~ 1 j) e- c( - s - r) dr 



+ / CT(-Y)((M - e - c ( s -*)) y -M Y + (G + e- c(s -^) Y -G Y )dr 
Jo 



First remark that when the underlying process F is a martingale, a = 1: indeed, from PIDE ( 3.4 ) , 
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Option H 


Moneyness 


Option value 


efficiency (8 true ) 


efficiency (6 mart ) 


Variance 
Reduction 


Call 


1 


4.90 


4.73 


5.03 


-6.34% 


Put 


1 


4.65 


4.70 


4.99 


-6.17% 


Call 


1.1 


1. 


6.02 


6.32 


-4.98% 


Put 


1.1 


8.63 


6.03 


6.31 


-4.64% 


Call 


0.9 


10.52 


3.58 


3.81 


-6.42% 


Put 


0.9 


2.33 


3.48 


3.69 


-6.03% 



Table 3: Pricing and standard deviation of hedged portfolio in the CGMY case 



we have 



da da 
dt ^ dz 
a(T, z) = 1 



(a(t,z + -/)-a(t,z)-j^(t,z)]v(dy)- inf \2irQa{z) + 7r 2 £a} 
V oz ) M<n 



On the other hand, from Definition 3.1 we have 



Qa(z) := / (e 



- 1) (a(t, z + 7(t, z, y)) - a(t, z)) v{dy) 



since fl, given in (2.2 1, is equal to zero (it is the drift of the process F which is now a martingale). 
From this, it is straightforward to deduce that the function a = 1 is the unique solution of PIDE 



(3.4 1. So that, when F is a martingale, one only needs to compute the function b. 



We now evaluate the loss of efficiency when using the martingale hedging strategies compared 
to the quadratic hedging strategies under the true historical measure. Our efficiency comparison 
criterion is the following: if H is the (put) option and Q true ^Q mart are; respectively, the optimal 
quadratic hedging strategy and the martingale strategy, then the efficiency is measured in terms 
of the standard deviation of the hedged portfolios: 



efficiency^'" 16 )' := Var H{F d>T ,t) 



0l™ e dF 



d,T,r 



(5.7) 



where x is the true optimal price given in (3.9). Similarly 



efficiency ( 



qmart\2 



:= Var (^H(F d ^ t ) ~ x mart - £ 



■d,T,r 



(5.8) 



where x mart is the price given in (3.9) when one uses the functions a and b computed in the 
martingale model., i.e. x mart is the risk neutral price of H. The variances are computed by Monte 



Carlo over 100000 paths using the rejection method algorithm described in Madan and Yor (2005 1, 



with 800 rebalancing dates in each path. The number of Monte Carlo trajectories used is limited 
due to the cost of the simulation algorithm. The trajectories of F d ,T,t are simulated using the 
true model in both cases. Table [3] summarizes the results of simulations, for t = 0. The numerical 
experiment proves that one loses efficiency when using the martingale hedging strategy. This is 
consistent with the fact that Q true achieves the minimum in problem (5.5) and outperforms the 
strategy mart . 



5.2 Numerical example: the Normal Inverse Gaussian process 



In this last paragraph we study the problem (5.5) when L in (5.1) is a Normal Inverse Gaussian 
process with parameters a,/3,8,u: L t ~ NIG(a, (3,5t,ui). 

Remark 5.2. The parameter a should not be mistaken for the parameter in Lemma \5.1\ We use 
this notation because it is standard in the literature. 
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We can write then 



05 



P 2 



t + 



yj{dyds) 



where J is a compensated Poisson random measure with intensity 

v(dy) = v{y)dy, v{y) = —r-K 1 (a\y\)e l3v , 



where K\ is the modified Bessel function of the second kind (Section 4.4.3 in Cont and Tankov 

1 r -(a-0)y v ,yj V^~°° 1 n -(oc+P)\y\ 



(2004 1). The Levy density v{y) satisfies 



v {y) 



*\y\ 



v {y) 



y—¥-\-oo 



\y 



3/2 



\y 



3/2 



Remark 5.3. The NIG is a infinite variation Levy process with stable-like behavior of small 
jumps, and since the Blumenthal-Getoor index is equal to 1, we cannot formally apply Lemma \5.1\ 
and Theorems \ 3. S\ [?Q| It is nevertheless a case of interest because the NIG model is popular among 
practitioners, and we shall see in the sequel that our numerical schemes yield good results for this 
model. 



We wan t to solve numerically the equations (|3.4|) and (|3.6|) where the model for Z is given in 
Lemma 



5.1 



and the source of randomness L in (5.1 1 is given by the NIG process introduced above. 
We apply the scheme (4.9| for the function a with maturity T = 1. Once again, the coefficients 

06 



H(t,z) =$ ($- L (z)) lu + 
j(t,z,y) =$ (^- 1 (z)+ye ct ) 



or 

— z 



2 







]e ct 4 


/[ 




Jr l 



$ ($ _1 (Y) + ye ct ) -z- ye ct $ ($ _1 (z)) v{dy) 



and the function D(t, z) introduced in (4.5 1 are computed by numerical integration over the points 
yi such that j(t, z, yi(t, z)) = iAz, or equivalently 



yi(t,z) := e 

By expanding 7 around zero we obtain 



(^- 1 (z + iAz)-^- 1 (z)) 



D(t,z) := 



y- K -i/2(t,z) 



7 (i,z,y)^(^)^e 2ct ($ ($- 1 (z)))^ 



V»+i/a(t,2) 



y 2 v(dy) 



U- X (z + (« + i)Az) - *-\z - (« + 



7T 



since, around zero, we have y 2 v(dy) ~ f + ^-y + 0(y 2 ). (See for example |Raible| dgOOOh ). We 
proceed as for the CGMY case, even if we do not have a priori results on the existence of a smooth 
solution. 

We consider again problem (5.5) for European options /, with maturity one week and delivery 
for the 7 days of the week as in (5.6). The parameters of the NIG process are u = 0.08, a — 6.23, 
= 0.06, S = 0.1027. The mean reverting coefficient c is taken equal to 0.19. 

We first analyze the convergence of the algorithm as function of the number of time steps and 
the mesh size N. Table [5] plots the approximated value of the hedging strategy for a call option 
with moneyness values 1 and 1.5 for different values of N and Nt- Compared with the CGMY 
case, the convergence as function of the space mesh size is somewhat faster, which is consistent 
with the theory (the NIG process has the Blumenthal-Getoor index equal to 1). 

For the remaining numerical experiments we suppose that M = 200, N = 400, Az — 2 x 10~ 2 
and the number of time steps equal to 800. The value k for defining the diffusion zone in the PIDE 
is equal to 3. 

We apply scheme (4.9) to obtain a numerical estimate of the function a (Figure 13 and the 
optimal control tt* (Figure [5]). Figure [6] shows the optimal price b for an at-the-money call option 
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Figure 4: The value function (t, z) H> a(t, z) for the NIG process. 



Pi curve 




Figure 5: The optimal control (t, z) M- 7r*(t, z) for the NIG process. 



and Figure [7] shows the function b for an at-the-money put option. The behavior of the different 
functions is roughly similar to that observed in the CGMY case. 

As in the CGMY case, we estimate the loss in the efficiency of the hedge when using a martingale 



model in terms of the standard deviation of the hedged portfolios as in (5.7)-(5.8). For F to be a 
martingale, we should have 



F, 



1 



d,T,t 



T+d 



V>(0, s) exp (M(s, t) + l(s)A t ) ds 
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B curve for call 




Figure 6: The value function b(t, z) for an at-the-money call option. 



B curve for put 




Figure 7: The value function b(t, z) for an at-the-money put option. 

where 



M(t,s) = - 




We already know that in this case we only need to compute the function b (a = 1 in the martingale 
case). Table [I] summarizes the results of simulations, for t = 0. The numerical experiment proves 
that using the martingale hedging strategy is inaccurate. The loss of efficiency is of order -20% for 
options with moneyness equal to 1. The magnitude of the non martingale model efficiency raises 
with the moneyness and can reach 100% for the Put option with moneyness equal to 1.5. 
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Option H 


Moneyness 


Option value 


efficiency^*™ 6 ) 


efficicncy((9 mQr *) 


Variance 
Reduction 


Call 


1 


4.199 


1.085 


1.316 


-21,29% 


Put 


1 


4.213 


1.087 


1.315 


-20,98% 


Call 


1.5 


0.120 


0.168 


0.212 


-26,19% 


Put 


1.5 


38.80 


0.175 


0.34 


-94,29% 



Table 4: Pricing and standard deviation of hedged portfolio in the NIG case 



Time steps (Nt) 


Space steps (TV) 


Price (moneyness 1.5) 


Price (moneyness 1) 


800 


800 


0.1201 


4.199 


400 


800 


0.1197 


4.274 


200 


800 


0.1188 


4.424 


100 


800 


0.1167 


4.718 


800 


400 


0.1215 


4.239 


800 


200 


0.1282 


4.505 


800 


100 


0.2028 


5.389 



Table 5: Convergence of the optimal price of the quadratic hedging strategy as function of time / 
space discretization parameters. 
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A Convergence analysis 

Proof of Theorem \4-2\ Under the additional assumptions of this theorem, 

D(t,z)= [ 7 2 (i,z,y)i/(dy) > rrn / y 2 v(dy)>m 1 { y 2 u(dy) ^ C(nAzf- a , Az -> 0. 



7r> ;o 71 i--f«+l/2)Ai 

where, throughout this appendix, C denotes a finite constant, which does not depend on any 
truncation / discretization parameters and whose exact definition may change from line to line. 
On the other hand, 

<Jh + Az l*M*>*)- 

k<\i\<M 

Using the mean value theorem and our assumptions on 7, we can show that 

IA(*. z )\< P- + C [ \v\v{dy) < fi + C\kAz\ x -<*. 

Therefore, one may choose a constant c which does not depend on Az (for Az small enough), such 
that for k > c, the weights a and /3 are both positive and given by (4.3)— (4.4). Throughout this 
proof we shall assume that such a choice of k has been made. 
Let 

/3 Z '*" = inf{i > n : Z*' tn + MAz > N Az or Z t z,t " - MAz < -NAz}. (A.l) 
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We have, 



a n (z 3 ) 



inf E 

7r 4 e[-n,n] 

tn, *3AN T 



n i+^(e Zt " ,j - 

i=n+l ^ 



1) 



1 



/ p -> i\ iV 

tn ' x 3 AN T 



inf E 

7r«e[-n,n] 



0*»-*i<N T 11 I 1 + 7r i( e 

II 1 + ^ -1) a(^,^') 



Let 



inf E 

iriG[-n,n] 



<-[jtn,*j >N , 



1)) a{t NT ,zl n N * 3 ) 

i—n+l 

Clearly, a n (zj) > a n (zj). On the other hand, by the Cauchy-Schwarz inequality, 



a„(zj) - a n (zj) < sup _ E 
7r;e[-n,n] 



i— n+1 



< ||g a |UF[/3*"'^ < 7V r ]i sup E 

7r t e[-n,n] 



1) 



The second factor can be bounded using the fact, that by Lemma A.l for every 7r e [—11,11]] 
E[(l + n(e 2 ^- z - l)) 4 ] = 1 + At [4w{»(t n , z) + J ( e 7(*-.*,w) _ i _ 7 (* n) z , 
+ 6tt 2 / ( e 7(W,s/) _ lf v (dz) + At: 3 [ ( e 7(*n,z,j/) - lfv{dz) + tt 4 / (e^*-"**) - lfv{dz) + £ 



where £ is the error term. By our assumptions, the terms in square brackets are bounded, and 
therefore, by applying iterated conditional expectations, 



sup E 

7r 4 e[-n,n] 



1) 



< (1 + CM) NT - n < e 



N T -n ^ „CATN T _ „CT 



(A.2) 



Together with the estimate of Lemma [A. 2| for the first factor, this yields 



\a n (zj) - a n ( Zj )\ < 5m}a?- flT ' ^Az > MAz + pT. (A.3) 
The remaining error is decomposed as follows: 



a n (zj) - a(t n ,Zj) 

E f inf 

7Ti£[-n,n],i= 



E 



=n+l,...,fc 



i— n+1 



inf E 

7r < e[-n,II],»=rj+l,...,A:-l 



fc-1 



i— n+1 



_i II (l + TiCe^""-^ -1)) a^.!,^/) } 



(A.4) 



2G 



The first term inside the brackets satisfies 



inf E 

7ri6[-n,n],-i=n+l,...,fe 



i=n-\-l 



W -1)) a(t k X7 Z] ) 



inf E 

7rie[-n,n],i=n+i,...,fe-i 



inf E 

7r fc e[-n,n] 



k-l 



z=n+l 



1) 



l + 7 r fc (eV'--V'-:L)) a(i fc , 



tk — l 



By Lemma A.l for every 7r S [—11,11], 



E 



l + 7r(e Z *»' -1)) o(t fc ,2£'") 



E 



l + 7r(e '* -1) o(t fc _i,Z'"^) 



*fc-i 



E 



da 



£t z 3 = 
tfc — 1 



dt 



where 



rtk da 

a(tk-i,z) + At£ w a(tk-i,z) + At / £ w — (i, z)dt + a(t k , z) - a(t k -i,z) + At£ k , 



\£ k \ < CAz' j - Q (^- Q + k ) +C (1 



\+T(y)+T 2 (y))u(dy), 

J\v\> M ,^ z 

Halloo 

and the operator is defined by 

C*f(t, z) = z)f'(z) + J{f{z + 7 (t, z, y)) - /(z) - /'(z) 7 (t, «, y))f (dy) 

+ 2tt ( M (t, z)/(z) + / (( e ^^) - i)/(z + 7 (t, z, y)) - /(z) 7 (t, z, y)>(dy) 
+ 7T 2 / ( e ^^>-l) 2 /(z + 7 (M,y)M^). 



From the regularity of a and the integrability conditions on the Levy measure, we deduce that 

da | 
dt 1 

2 



C-x%{t, z) is uniformly bounded on it <E [—11, II], which means that 



E 



l+ir(e z *> ~ z **-i -1)) a(t fc ,^-'*) 



<9a 



with 



|£ fe |<C£, £ := At + Az^- Q (^- Q + + 



a(t k -x,z)+Ab \ C n a(t k ^i,z) + —(tk-i,z) \+At£ k 



{l + \y\+T{y)+T 2 {y))v(dy), (A.5) 



l«l>i 



Using the equation (3.4) satisfied by a, we finally get 



inf E 
w6[-n,n] 



a(t k -\,z) 



< CAt£. 
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Plugging this estimate back into (A.4) yields 



a n (zj) - a(t n ,Zj) 

Ni 



< T ( inf 

_ y 7rie[-n,n],i=n+l,...,fc-l 



E 



fc=n+l 

inf E 

7riG[-n,II],j=n+l,...,fe-l 



k-1 



k-1 



i n 



i)i {a(t k - u z£*) + ctee) 



i— n+1 



~-t T , ,2.- -~-t„ 



^*-»i>jk-i II fl + ^(e Zti ' 3 ^ -1)) a(i*_i,^) } 



which implies: 

a n (zj) - a(t n ,Zj) < CAtE ^ 



sup E 

k=n+l T«e[-n,fi],*=n+l k-1 



k-1 



2=71+1 



n 



1) 



The expectation can be estimated as in (A.2), and we finally get 

a n (zj) - a(t n , Zj) < C£. 
The upper bound can be obtained in a similar manner. 
Proof of Corollary 



□ 



4-6 We first estimate the error of approximating the numerator and the de- 



nominator in 1 4.11 1. For the numerator, we get: 
E 



e '■n.+ l 



* -1 ) a n+1 {ZlZ>*l) 



AtQ t a(t n , z 3 



< 



+ 



e 3 - 1 (««+! ) - a(t n+1 , 



e f " +1 J - 1 (a(t n+1 , Z™1) - a(t n , Z£%)) 



AtQ t a(t n ,Zj) 



+1 ''-I U(<„,^) 



Using Lemma A.l we can show that the second term is bounded by CAt and the third term is 
bounded by CAtE, with £ defined in (A. 5 1. The first term makes the main contribution to the 
error, which can be estimated using Theorem |4 . 2 1 and the Cauchy-Schwarz inequality: 



E 



* - 1 ) (a n+1 (zi:^) - <*(*»+!, 



< CSE 



z'" 3 



< CVAtS 



Similarly, for the denominator we have the estimate 



E 



J -l a, l+1 (z;;;7) 



- AtQ t a(t n , Zj) 



< CAtE. 



Using the fact that for b, b' > and all a, a' , 



a a 
V~b 



< 



\a! - aj KH^-bl 
6 + W ' 



we obtain the estimate 



E 



^+^-1)^+1 (<";?) 



Q t a(t ni Zj) 



E 



< 



Cy/AT£ CAti£ 



Atg t a(t n ,Zj) 



Atg t a{t n , Zj ) At 2 g t a(t n , z 3 ){g t a{t n , Z] ) - C£) 



< CAt~*£ 
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since Q t a is bounded from below (this follows from Theorem 2.3 1. We conclude by observing 
that projecting both the optimal strategy and its approximation on the interval [—11, II] does not 
increase the error. □ 



Proof of Theorem 4-8 As in the proof of Theorem 4.2 we may and will assume that a and /3 are 



positive and give n by (4.3 1— (4.4). To deal with the domain truncation, we also, similarly to the 



proof of Theorem 4.2 take /3 2 >*™ as in (A.l), and use the representation 



b n ( Zj )=E 



1/3'™ < z i <N T 



n 

i— n+1 



A'T 



~*-p tn ' z j>N T 



i — n+1 



1 + 7r*«(Z t 7"')(e * - 1) q "(t^ JANT ,Z^2 ) 



JJ ( 1 + 7r*«(^^)(e^ B ' V -^ ~ 1) ) f(Zfc?) 



b-{ Zj ) = l^ j>NT _ x J] [l+^iZtmie ^ -1) )f{Zllt i ) 



Letting 



the remainder can be estimated as follows: 

\b n {z ) - b n ( Zj )\ < \\q b \\ocnP tn ' Zi < N T ]*E 



where the last inequality uses the fact that tt* is the minimizer of the expectation in the first line, 
and that the expectation equals one when substituting the value ir* = 0. Therefore, \b n (zj) — b n (zj)\ 
admits the same bound as \a n (zj) — a n (zj)\, given by (A.3|. 
The remaining error is decomposed as follows: 



N T 



b n (z j )-b(t n ,z j )= ]T 



k=n+l 



k-l 



i— n+1 



x E 



l+n *k { zi : *)( e z *i' i - z *i-i - 1) ) b(t kl zi7 Z] ) b{t k ^ziz z :) 



Similarly to the proof of Theorem |4.2[ by Lemma |A.1| we deduce that 



E 



l + 7r* k {z)(e <<= 



i))b{t k ,zl;^)-b{t k - u zlii z i) 



db 



At (B t - A t )b{t k - U z) + n**(z)Q t b(t k - U z) + —(t k - u z) > + At£ k 



dt 



= At(ir* k (z) - n*(t k - 1 ,z))Q t b(t k - 1 ,z) + At£ k 



with £ k satisfying (A.5), where in the last equality we used equation (3.6). We conclude by sub- 
stituting the estimate for \ir* k (z) — ir*(t k _i, z)\ given in corollary (4.6), and using the fact that 
Qt&(i/t-i, z) is bounded due to the regularity of b. □ 

Auxiliary lemmas 
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Lemma A.l. Let f be 4 times continuously differentiable with bounded derivatives. Then, 



E[f(Z?±)]-f(z) 



At 



V(t n ,z)f'(z) - / (f(z + j(t n ,z,y)) - f(z) - f'(z)-t(t n ,z,y))v{dy) 



< CAz 3 ~ a (k 3 - a + k 1 - ) + C 
E[(e 2 ^- z -l)f(Z^ :i )\ 



(l + \y\)u(dy) 



(A.6) 



' llTyllc 



At 



- fi(t n , z)f(z) - / (( e 7(*».*,») - l)f( z + 7 (t n , z, y)) - f(zh(t n ,z, y))v{dy) 



< CAz 3 ~ a (k 3 - a + k 1 -") + C [ 

J\v 



r(y)v(dy) 



(A.7) 



' ll-Yj/llc 



_ , (e7(t „, z , a) _ l)Pf{z + 7(tnj Z) y))i/(dy) 



< C*Az 3 - Q (fc 3 ~ Q + /c 1 -") + C / 



T p {y)u{dy), p = 2,3,4, 



(A.8) 



■ Hylic 



where t is defined in Assumption \2. 



Proof. We begin with (A.6). By definition of Z, the expression inside the absolute value on the 
left-hand side may be decomposed into the following terms: 

/(* + Az) + f(z - Az) - 2f{z) - /" (z)Az 2 



2Az 2 



-D(t n ,z) 



(A.9) 



f(z + j(t n ,z, y)) - f(z) - f'(z)j(t n , z, y) - -f"(z)^(t n , z, y) \ v{dy) (A.10) 



f(z + Az)-f(z-Az)-2f'(z)Az 
2Az 



n 3 (t n ,z) 



(/(« + 7(*n, z, y)) - /0) - 7(*nj ^, y)f{z))v{dy) 



E 



i:i/i(t»,«)6fla ^'-VaC'*) 
rVi+i/2{t,z) 
+ E / 

where we denote 



(F(z, 7 (t n , z, j/j)) - F(z, 7 (i n , z, y)))v{dy) 
(F(z, 7(t„, z, t/i)) - F(z, 7 (i„, z, y)))y 2 v(dy), 



(A.ll) 
(A.12) 

(A.13) 
(A.14) 



^(z,y) = /(z + y)-/(z)-y/'(z) 

and /=:. 

7 (t,z,y) z 

The function F clearly satisfies 

H^Hoo < 2||/'||oc, ll^'lloo < liriloo. 

As for the function F, we can write it as F ~ u(z, y)v 2 (t 1 z, y) with 

f{z + y)-f{z)-yf{z) 



u{z,y) 



f"(z + 9y)(l 



and v(t,z,y) = —— - = / 7y (6> 7 ^y))^. 

7 \ T , z ,y) Jo 
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It is easy to see that under our assumptions the functions u and v are bounded together with their 
first and second derivatives. It follows that ||i^||oo < oo and ||oo < oo. 
Remark also that by the mean value theorem and our assumptions on 7, 



^2 Ui(t,z)j(t,z,yi(t,z)) 



<C I \-r(t,z,vMdy) 

fiiUfi 2 



for some constant C which does not depend on truncation / discretization parameters. 
The terms in ( |A.9|A.14 1 admit the following bounds: 



|(A.9l| < 



ll/(4)|U -Az 2 / ^t n ,z 7 yMdy) <M^ Az 2 hy \\l 



24 



24 



(Hl/2)Ai 



y 2 v(dy) 



< CAz 2 (kAz 



,2-Q 



|(|Al0l)| < 



ll/ (3) ll 
6 

ll/ (3) l 



^(t^z^Hdy) < C(kAz) 3 



l dAH| )| < MJaz 2 (/1+ / 



\j(t n) z,y)\u(dy)) 



<M^ Az2(A + n 7y || f 

J\v\> [ - 



\yHdy)) < CAz 2 (kA Z y~ a , 



WfyWa 



|(|Al2|| < 2||/|| 00 i/(» 3 ) + ll/lloo / \j(t,z,yMdy)<C 

r(»+l/2)Az 



{l + \y\)v(dy), 



\y\> 



|(A.13)| 



J2 

i:y i (t n , z )en 2 J ( i - 1 / 2 ) Az 



(F(z, - F(z, iAz)) V (Z.-$l ^ 



< E 

i:yi(t n ,z)£Sl2 



Az 2 



\F„. 



Vi+1/2 



7v(7" 1 (0) 

(j+l/2)Az 



Vi-1/2 



<CAz 2 [ v{dy) + II^IU' 

.... „ .N-r, 



-A x ^(7~ 1 (0) , r 

(i-l/2)A* 7jA7 HO) 



d( 



r{i + l/2)Az 


d K7" 1 (0) 


/(i-l/2)Az 


rfC7y(7" 1 (0) 



< CAz" 



i/(dj/) 



n 2 



^ k'(y)|d2/) < CAz 2 



Similarly 



|(A.14)| 



E 



(i+l/2)A* _ „ ufa-ifO) 

(F(z, - F(z, »Az)) 7 ' 1 (C) 2 ; V =j£fe d( 



< 



l :y 1 (t„, 2 )G0 1 J ( l - 1 /2)A Z 

Az 2 



E 



F 



Vi + l/2 



y 2 u(dy) + \\F v \ 



Vi-l/2 



7y(7" 1 (0) 

(i+l/2)Az 
(•i-l/2)Az 



d_ 



7" 1 (C) 2 



K7-H0) 
7,(7- 1 (0) 



< CAz 2 



+ / |vK<i|/) + / |j/V(j/)|dy <CAz z (kAz) 



"1 



We finish the proof of (A.6) by observing that given that kAz is small (this is the small jump 
truncation level), the leading contribution is made by terms (A. 10 1, (A. 11 1, (A. 12 1 and (A. 14 1. 
We next prove (A. 7). The expression inside the absolute value on the left-hand side may be 
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decomposed into the following terms: 

(e Az - l)/(z + Az) + (e- Az - l)f(z - Az) - (f(z) + 2/'(z))Az 2 

2Az 2 



D(t n ,z) 



(A.15) 



(e Az - l)f(z + Az) - {e- Az - l)f(z - Az) - 2f(z)Az , 



f2 3 (*n,z) 



(e 7(M,y) _ + 7(tn) Zj y)) _ f(zh(t n , z, y) - l -{f{z) + 2/'(z)) 7 2 (t„, z, y)| 

(A.16) 
(A.17) 
(A.18) 



2Az 



fi(t n ,z) 



(( e 7(W,lO _ + 7 ( tn; Zj y )) _ 7 ( 4nj Zj y)f( Z ))v{dy) 



E 



J/i+l/2(t.«) 



(in, z)efi 2 ,/ »<-i/2(*.' s ) 
rVi+i/a(t,z) 

+ / , , 

where now the functions 



(G(z, 7(t n , z, j/j)) - G(z, 7(< n , z, y)))v(dy) 
(G(z, 7 (i„, z, j/i)) - G(z, 7(f„, z, y)))y 2 v{dy), 



(A.19) 
(A.20) 



G(z,y):=(e»-l)/(« + y)-i//(z) and G(z,y):= 



(e»-l)/(z + y)-y/(z) 

7 _1 (M,i/) 2 



satisfy 
and 



|G;| + |G"|| <C(e« + l), ye 



|G;| + IG'JI < G, |y|<l. 
The latter estimate, in particular, follows from the Taylor formula representation 
(ev -l)f(z + y)-yf(z) 



ir 



{e 9 y(f(z + e y )) + 2e y(f(z + e y )) + (e»» - l)(f"(z + 6y))}(l - 



The different terms can be estimated in a manner, similar to the first part of the proof: 
|flAT5l)l < GAz 2 (fcAz) 2 " Q , 



|(|Al6l)| <G(fcAz) 3 -", 
T(y)v(dy), 



|(A.17)| < GAz 2 (/cAz) J 



\ (aM \ <c [ 

|( [A~19| | < GAz 2 ( / ( e 7(W,2/) + l)v{dy) + f (e^"'^ + l)\i/(y)\dy) < GAz 2 
\Jn 2 Jn 2 / 



|([A~20|| < CAz 2 (kAz 



,1-a 



Finally, (A.8) can be proven in a similar manner. For example, for p = 2, we decompose the 
expression inside the absolute value on the left-hand side into the following terms: 

(e Az - l) 2 /(z + Az) + (e- Az - l) 2 /(z - Az) - 2/(z)Az 2 

2Az 2 V[t n ,z) 

( {e y(t,z,y) _ 1)2/(z + l{fni ^ y)) _ /(z ) 7 a (tnj Zt y) \ v[dy) 

n Q (*„,*) L J 
(e Az - l) 2 /(z + Az) - (e- Az - l) 2 /(z - Az) , 



2A2 



-fi(t n ,z) 



(e^^-l) 2 f(z + 7 (t n ,z,y))u(dy) 



rVi+-L/2(t,z) 

■■ yi {t n ,z)en 2 Jvi-i/z&z) 
rVi+i/2(t,z) 

:: yi (t n ,z)eni ^-i/2(*' z ) 



(#(z, 7(i„, z, - i?(z, 7(f„, z, y)))u(dy) 
(H(z, j(t n , z, y»)) - AT(z, 7(t„, z, y)))y 2 v(dy), 
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with 



H(z,y):=( e y-l) 2 f(z + y) and H(z, y) (< '" 



7 1 (t,z,y) 2 

□ 

Lemma A. 2. There exists a constant C < oo which does not depend on the truncation / dis- 
cretization parameters such that for all n £N and all K > fmAt, the approximating Markov chain 
Z satisfies 

™r ,/„ i z 2 + CnAt 

P max \Zf! to \ > K < — — t-^-. 

L o<i<™ 1 *« 1 - J - (if - £nAi) 2 

Proof. By construction, the approximating chain Z satisfies 

E^J-ZtJ = //(**, £ t JAf 

Var[Z ti+1 |Z f J =^(t l ,Z f JAt-/i 2 (t 4 ,Z t JAt 2 + At ^ u^uX^iuX^A))- 

i:i/i(ti,z ti )eniuna 

Therefore, the process M ti = Z|? °— X^o<j<i Ztf°)At i s a discrete-time martingale. By Doob's 
martingale inequality we then get 

E\M 2 ] 

Pf max l^'* 1 > itl < Pf max |M t . I > K - uAtn] < - 1 *" J NO 

l o<i<n l u 1 - J - L o<i<n' p J - (K - fiAtn) 2 

^+E[Eo< 3 <nVar[% t ^f"]] 
(if - jlAtn) 2 

< ^ 2 +EEo< J -<»Pfe^ tj ) + E ;: , ig o 1 uo 2 7 2 fe,^,wMfe,^-))]At 

(if - //Ain) 2 

We finish the proof by applying the mean value theorem and using the boundedness of the deriva- 
tives of 7. □ 



B Proof of Proposition 5.1 



Proof. Before we start, remark that the function A H> F<i,t{A) is strictly increasing, so invertible, 
and infinitely diffcrcntiablc: in particular 



J^ +d ^(0,s)l(s)e l ^ A ds 



<$>"(A) 



/; + >(0,s)e'( s )^s 

f^ +d i(>(0, s)l 2 (s)e l ^ A ds\ (/J +< V(0, s)e z(s)A ds) - (j™ ^(0, s)l(s)e l ^ A di 

^ +d ^(0,s)e'^ A d s y 
from which we deduce 

e -c(T+d) < < e -cT and e -2 c (T + d) _ e -2cT < < e -2cT _ e -2c(T+d) 

From Ito's formula, we obtain 

dZ t = (^>'(A t )e c \ + J {*{At- + e ct y) - $(A t _) - ye c *$'(A t _)) i/(d tf )) <tt 
+ y ($(A t _ + e ct y) - $(A t _)) J(dydi) 

or equivalently 

dZ t = //(i, Z t )dt + / 7(t, Z t _,y)j(dydt) 
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We can now prove that fi and 7 verify the Assumptions |2.l| We detail the computations only for 
the function 7, since similar computations can be done for /1. First we remark that z 1— > 7(4, z, y) 
is diffcrentiable and we can compute this derivative to obtain 

0, 7 (t, z,y) = -l + ($'($- 1 (z)))- 1 $'(<i>- 1 (z) + ye ct ) 

fi 



so that 



\dzl(t,z,y)\ = 



=e ct y($'($- 1 (z)))- 



e^ y ($'($-i(^)))- 



$"($- 1 (z) + re ct )dr 



$"($ _1 (z) + re c *)dr 



<Me cT (inf ^(A)!)- 1 ||$"|L < e cT e- c ( T+d ) ll^'IU |y| 



<|y|e cT e c ( T+d ) fe- 



-2cT -2c(T+d) 



< e cd \y\ 

From the bounds on the first and second derivative of $ we obtain sup t z \d z "j{t, 2, J/)| < e cd |?/|, 
which gives us the function p introduced in Assumptions 2.1 Again by the definition of $ in (5.3) 
we have 

exp( e - c(T+d) y) - 1 < e 1 ^'"^ - 1 < e v - 1 

if y > and the inverse inequality stands in force if y < which yield sup t z \e y ^ t,z ' y ' — 1| < |e y — 1|. 
According to the definition of the function r given in Assumptions |2.1| and the estimations above 
we deduce that 



t (y) := max fsup ^7 (t, z, y) |, 



p(y) 



If follows then that Assumptions 2.1 [C, I, L] hold true. For Assumption 2.1 [ND] we have, from 
the definition of 7 

(erf*.*.»> - I)' > (exp(e- c( - T+d 1y) - lY 
so then, for some positive M > we have 

T(y) := { inf (e^*'*'^ - l) i/(dy) > / inf (exp( e - c(T+ti) y) - l) v{dy) 
Jb M \ / Jr M \ / 

lyl 1 -^^)^ > 



>M 



\v\<e 



since g(0 + ) and g(0 ) are strictly positive, we can select e small enough and obtain 



T(y) >M f 

J\v\<* 



t/l 1 "^ > 



We can derive 7 w.r.t y to obtain 



7„(M,!/) =e ct &{$-\z) + e ct y) 
7 w (M,y) =e M $"(r 1 (z) + e d y) 



2.1 



RGj] holds true. For Assumption 
an upper bo und to prove that indeed z 
Assumption 



2.1 



so then e r c ( T + d ) < \j y (t, z,y)\ < e cT and \j yy (t,z,y)\ < e 2cT , which proves that Assumptions 

[RGjii], one can differentiate j y w.r.t. z and give for it 
J y (t, z, y) is Lipschitz continuous uniformly in t, y. The 
[RGii] does not hold true since trivially "f y (t, z, y) — e c *<f>'($ _1 (z) + e ct y) 7^ 1 . 



2.1 



The last thing we need to prove is the condition 4 assumed in Theorem 4.2 i.e. that 7 is 3 times 
differentiable w.r.t. y with bounded derivatives. From the definition of 7, this is equivalent to prove 
that $ is 3 times differentiable with bounded derivatives. Let us introduce 



Pi(A) 



T+d 



ip(0, s)l l (s)e l ^ A ds 
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so that = p 1 (A)/p (A) and = (p 2 (^)Po(^) - Pi(^) 2 )/Po(^) 2 - By remarking that 

= pi + i(A), we can differentiate $" to obtain 

^ui) {A) = P3(A) P0 (A) +P2{A)pM) 2 Pl (A)P2(A) 2 #, {a) # {a) 

with $' and $" bounded as already proved. With the same type of computation it is straightforward 
to prove that $( W )(A) is also bounded. 

□ 
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